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I. INTRODUCTION 



One of the main goals of modern nuclear physics is to study the behavior of nuclear 
matter under extreme conditions. At low temperatures and densities the quarks and gluons 
as the basic constituents of strongly interacting matter form hadrons due to the confinement 
mechanism. In addition, the appearance of rather light mesons (pions and kaons) signals the 
existence of a spontaneously broken symmetry, the chiral symmetry. In fact this symmetry 
is approximately realized in the QCD Lagrangian. Another important hint that chiral 
symmetry is spontaneously broken in the vacuum state is the absence of chiral partners 
with equal masses. In a chirally symmetric state chiral partners would have the same 
mass. This concerns for example the isovector vector meson p and its much heavier partner, 
the isovector axial- vector meson a±. It is expected that at high enough temperatures and 
densities confinement is lifted and chiral symmetry restored. High energy heavy-ion collisions 
are dedicated to the creation of this new state of matter, the quark-gluon plasma (QGP) 
|Ij]. Unfortunately, even if such an ultra- hot system of quarks and gluons is created only 
its decay products — which of course are hadrons and not deconfined quarks and gluons - 
can reach the detectors. Thus, the proof for the existence of this new stage of matter has to 
be performed rather indirectly. Besides the observable hadrons also photons and dileptons 
are radiated from the hot fireball. These particles deserve special attention since they do 
not suffer from strong final state interactions. Therefore, once created in the high density 
region they are capable to carry information from that region to the detectors. Altogether, 
the challenge is to find unambiguous signs that (part of) the observed spectra of hadrons, 
photons, and dileptons are caused by the transient existence of the QGP. Clearly, to prove 
the existence of the QGP it is necessary to show that the spectra cannot be explained by a 
hot fireball made out of conventional interacting hadrons. This task is especially complicated 
by the fact that there is no straightforward derivation of hadronic Lagrangians from QCD as 
the underlying theory of strong interactions. Therefore it is not a priori clear how far one can 
trust in-medium calculations with hadronic Lagrangians as their parameters are adjusted to 
the description of vacuum processes. Connections between hadronic models and concepts on 
the one hand side and QCD or QCD based models on the other hand side are therefore very 
welcome. The QCD sum rule method provides such a link. Before we sketch its basic 
concepts, however, we want to dwell for a moment on the in-medium properties of hadrons. 

For temperatures and densities below but near to the critical values which mark the 
transition to the QGP it is plausible to expect that already there the properties of the 
involved hadrons like e.g. their masses and decay widths get modified. Especially the aspect 
of chiral symmetry restoration is interesting here. The properties of chiral partners should 
start to approach each other and finally become identical in the chirally symmetric phase. 
Concerning p and a\ mesons possible scenarios are e.g. discussed in |4]-|| for the case of finite 
temperature. In principle one can distinguish three types of possible phenomena (which do 
not exclude each other): 

a) Mass shifts: The masses of p and a\ might approach each other. One has to distinguish 
in which way this actually happens: The masses might meet at a value somewhere in 
between their vacuum masses (and possibly drop together afterwards). It is however 
also possible that the masses of both mesons drop and finally (approximately) vanish 
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at the point of chiral symmetry restoration. 



b) Peak broadening: From the experimental point of view the p (a\) meson shows up as 
a peak in the vector (axial-vector) channel. In a medium the peaks might get broader 
(maybe without a change of the respective peak positions, i.e. the nominal masses) 
until the melted spectra in both channels become degenerate. 

c) Mixing: The distinct peaks might maintain (maybe without shifts or broadening), but 
the a\ peak shows up with increasing height in the vector channel and vice versa. 

In any case, the spectra in the vector and axial-vector channel become degenerate when 
chiral symmetry becomes restored. 

In fact, the p meson is supposed to be a good candidate to search for a sign of chiral 
symmetry restoration. The reason is that it has the quantum numbers of the photon. 
Therefore the p meson can decay via a virtual photon into a dilepton pair. If this decay 
happens within a hot and dense medium the dileptons contain information about the in- 
medium properties of the p meson. Therefore, in principle the possible scenarios discussed 
above or a mixture of them should leave their marks in the dilepton spectra. Indeed, the 
HELIOS and CERES collaborations have reported medium modifications in the dilepton 
spectra in the invariant mass range around the p meson mass 0,0. Whether the observed 
spectra can be explained within a conventional hadronic scenario |§ or whether one has 
to include medium modifications induced by chiral symmetry restoration []10l,|TT]l is still a 



matter of discussion (cf. also fll2j|l3l and references therein). Also for the study of possible 
in-medium changes of hadronic properties a closer connection between hadrons and QCD is 
desirable. 

The QCD sum rule approach has the merit to relate certain low energy quantities - 
which so far are not directly accessible by QCD — with high energy expressions which can 



be calculated by the operator product expansion H] in terms of quark and gluon degrees 
of freedom. Non-perturbative effects are encoded in the appearance of various quark and 
gluon condensates. In the following this method is applied to vector and axial-vector mesons 
placed in a cold medium with finite nuclear density. To clearly work out the modifications 
when changing from vacuum to a medium we also discuss the vacuum sum rule analysis for 
p and ai in some detail. Before sketching the basic ideas of the sum rule approach we review 
the present status of in-medium analyses for vector and axial-vector mesons especially in 
the light of the possible scenarios of in-medium changes discussed above: Concerning finite 
temperature T it has been shown that at 0(T 2 ) and neglecting the pion mass only mixing 
occurs. A systematic study beyond this linear (pion) density approximation is complicated 
by unknown non-scalar higher twist condensates |15|-[17| . For finite nucleon density previous 
analyses have restricted their attention to the vector channel. In the first analyses |l8|Jl9[| 
only a possible mass shift for the p meson has been taken into account, i.e. the possible 
scenarios of peak broadening and/or mixing as mentioned above have been excluded by 
hand. In this case it was found that the p meson mass would drop in a nuclear medium. 
However, it has been shown by the authors of |20| that their specific hadronic model also 
fulfills the sum rule (cf. also |f2T| , p2f ). This model predicts peak broadening for the p meson 
and basically no mass shift. Subsequently, a systematic study revealed that independently 
of the chosen hadronic model the sum rule for the p meson for finite density is in accordance 
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with a specific mass- width correlation ||23|| : For low width the mass has to decrease. If, 
however, the mass stays constant — or even rises — the width has to increase. The sum 
rule does not have enough predictive power to fix both the mass and the width of the 
vector meson. We will come back to that point below. To the best of our knowledge, the 
third possible in-medium modification, the mixing phenomena, has not yet been included 
in a systematic sum rule analysis for the vector axial-vector system for finite density. The 
purpose of the present work is to treat the properties of p and a\ on equal footing, allowing 
for mass-shifts, peak broadening, and mixing. 

Before we turn to the specific sum rule analysis for the chiral partners p and a\ we discuss 
some important aspects of the QCD sum rule approach focusing especially on in-medium 
situations. Recall that our basic motivation was to describe various spectra of heavy-ion 
collisions by hadronic models in absence of a QCD description based on first principles. 
Even when only hadronic models are capable to calculate observable quantities one can 
imagine that it is possible to find other quantities which can be reliably determined both 
within the hadronic framework and in terms of quark and gluon degrees of freedom. In this 
way one obtains predictions for hadronic parameters like masses and coupling constants or 
cross-check for hadronic models. Concerning the QCD sum rule method such quantities 
are specific correlators (see below) calculated in the deep space-like region, i.e. for large 
momenta q with q 2 <C 0. For very large Q 2 = —q 2 QCD perturbation theory becomes 
applicable. Proceeding to (somewhat) smaller values of Q 2 non-perturbative corrections 
appear. They can be expanded in a power series in 1/Q 2 , called the operator product 
expansion (OPE): 

Y — . (1.1) 

In the coefficients c n the famous quark and gluon condensates enter. One can imagine the 
series ( |1 . 1| ) as a separation of the hard (denominator) and soft (numerator) scales of the 



problem (cf. e.g. [24| and references therein). In the numerator the non-perturbative effects 
enter. In practice, only the first few coefficients in ( |1 . 1| ) can be determined. Of course, this 
does not matter as long as Q 2 is large enough. Thus, the crucial question is for which values 
of Q 2 one can trust the truncated series. If we want to learn something about a hadron with 
mass nth it turns out that Q 2 has to be of the order of m\. To get an order of magnitude 
estimate for the coefficients c n we have to ask about the typical scales for non-perturbative 
effects. Let us discuss step by step the different cases of vacuum, finite temperature and 
finite baryon density. For vacuum the typical scales are Aqcd and the current quark masses. 
The up and down quark masses have only a few MeV and are therefore negligibly small. The 
strange quark mass and Aqcd are between 100 and 200 MeV. On the other hand, the typical 
hadron masses are of the order of 1 GeV. Therefore, one might expect that the sum rule 
analysis leads to reasonable results. [] Of course, the masses of the much lighter pions and 
kaons cannot be determined. Unfortunately this optimistic picture is not completely true. 



1 We restrict our considerations here to hadrons made out of light quarks. The masses of the heavy 
quarks have to be regarded as part of the hard scale . 
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In fact, there might be non-perturbative effects which introduce an additional hard scale, 



like e.g. instantons p5fl . In this case the series (|1 . 1| ) would break down for the interesting 
values of Q 2 . It seems, however, that the influence of such effects on the p and a\ sum 
rules is not important. We therefore assume throughout this work that the OPE works for 
the vector and axial-vector channel. Nonetheless, this consideration shows that at present 
the QCD sum rule approach cannot be directly justified from QCD without any additional 
assumptions. Therefore it should merely be regarded as a QCD based model and not as 
QCD itself. Turning to the case of finite temperature involves new scales. It is common 
practice to approximate the lowQ temperature medium by a pion gas. Therefore the new 
scales are the temperature and the pion mass. Also these quantities are of the order of 
Aqcd- Of course, T might also be lower. Therefore the previous considerations apply also 
here. The case of finite baryon density pn is more complicated. Here one approximates 
the medium by a Fermi gas of nucleons. New in-medium scales are the Fermi momentum 
and the nucleon mass. While the former is reasonably small e.g. for saturation density of 
nuclear matter, the latter is of the order of 1 GeV. Since the nucleon mass enters the series 
( |1 . ip in the numerator it becomes questionable whether the OPE still works (cf. also the 



discussion in |26] and the successive comments |p7 ,p8fl ) . Full clarification of this question 



requires the determination of all coefficients c„ which would be equivalent to solving QCD 
in the non-perturbative low energy domain. This is of course out of reach. For our case 
at hand there is, however, a class of contributions to the OPE which can be determined to 



all orders, namely the twist-two condensates [lS,iy,El,p9]. In fact their contribution to the 



coefficients in the low density approximation is given by 

C ist ^° = a n m%- 3 p N , (1.2) 

i.e. powers of Q 2 in the denominator which in an optimal situation should suppress higher 
order contributions in ( |1 . 1| ) are compensated by powers of m 2 N in the numerator. Thus the 
class of twist-two contributions shows exactly the unpleasant feature discussed above. In 



L2|) the dimensionless quantities a n can be determined from the parton distributions in a 



nucleon |l8| , p!9|p^ , |29[| . Fortunately it turns out that a n is strongly decreasing with increasing 



n such that the higher dimensional contributions of the twist-two condensates can safely 
be neglected |24,29||. This is a hint that the OPE still works in the case of finite nuclear 



density. Of course this is not a proof for the validity of the OPE. Throughout this work we 
assume that the OPE works. In spite of these obvious problems inherent to the QCD sum 
rule approach for finite density we regard the analysis presented in the following as useful in 
view of the possibility to learn something about the in-medium properties of hadrons from 
an approach which deals with the fundamental degrees of freedom of QCD. Nonetheless we 
stress again that the QCD sum rule approach — especially for the case of finite nuclear 
density — is not as fundamental as QCD. 

The article is organized in the following way: In the next section we derive the in-medium 
Borel sum rules for p and a\ which we will use throughout this work for any quantitative 



2 At high temperatures it is not reasonable to deal with hadronic degrees of freedom. If one wants 
to learn something about hadrons low temperature expansions are appropriate. 
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statements. In Sec. [Ill] we make a detour to discuss a different type of sum rule, namely 
the finite energy sum rule. This will yield a qualitative picture what one has to expect from 
a sum rule analysis and how much information one can get. In Sec. |IV] we introduce our 
hadronic parametrizations which are used to analyze the sum rules. Results for p and a\ 
in vacuum are presented in Sec. [V|. These results serve as a reference frame with which 
we can compare the succeeding in-medium results. Sec. |VT] is devoted to the discussion of 
the mixing phenomena while the in-medium results are presented in Sec. |V11| . Finally we 
summarize and discuss our results in Sec. [VIII . 



II. THE CURRENT-CURRENT CORRELATOR AND THE BOREL SUM RULE 

The relevant quantity to look at is the covariant time ordered current-current correlator 

IM?) = i Jd'xe^iTj^UO)) . (2.1) 
For the p meson channel is the isospin-1 part of the electromagnetic current, 

= o ( Ul ^ u ~ ■ (2-2) 

This current-current correlator enters e.g. the cross section of e + e~ —>■ hadrons (see below). 
For the a\ meson channel we have to deal with the corresponding axial-vector current 

it = o (^75« - dj^j 5 d) . (2.3) 

The expectation value in Q2.1| ) is taken with respect to the surrounding environment. 
We study here, first, vacuum and, second, an (isospin neutral) equilibrated homogeneous 
medium with finite nuclear density and vanishing temperature. In the medium Lorentz 
invariance is broken. All the formulae which we will present in the following refer to the 
Lorentz frame where the medium is at rest, i.e. where the spatial components of the baryonic 
current vanish. For simplicity we restrict our considerations to mesons which are at rest with 
respect to the medium. For the vacuum case we can choose the rest system of the (axial-) 
vector meson without any loss of generality. 

In the following the formulae without an explicit V or A index are valid for both vector 
and axial-vector channel. The correlator ( |2.1| ) has the following decomposition (valid for 
mesons at rest!) 

IW?) = %QuR(q 2 ) ~ 9,uIl isotT (q 2 ) ■ (2.4) 

In the following we concentrate on R(q 2 ). In the vector channel one has II lsotr (g 2 ) = q 2 R(q 2 ) 
since the current is conserved. We prefer the use of R instead of n isotr since it has been 
shown in [|T§| , |3"U| , |3"T| that the Borel sum rule (see below) is rather unstable for the latter 
quantity. The divergence of the axial-vector channel is solely determined by the pion decay. 
Hence we would not learn anything new about the a\ by studying n isotr in addition to R. 
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Concerning e.g. the dilepton production one is interested in the values of the dimension- 
less quantity Ry{q 2 ) in the time-like region q 2 > 0. The reason is that Ry is related to the 



cross section e + e — > hadrons with isospin 1 via [[32 

<r /=1 (e + e~ — > hadrons) 



a(e + e — ► yU + /i" 



127tIiilRv . (2.5) 



At least for low energies the time-like region is determined by hadronic degrees of freedom. In 
principle, there are two possibilities to describe the current-current correlator. First, guided 
by an educated guess one might use a simple parametrization with some free parameters. 



Second, one might use a hadronic model, e.g. for vector mesons ||13P0|^liP3| -p9[| using one 
or the other form of vector meson dominance. In the following we will explore the first 
possibility and figure out which constraints for these free parameters are provided by the 
QCD sum rule approach. For the a\ we proceed completely analogously. We denote the result 
for R in the time-like region by i? HAD . On the other hand, the current- current correlator 
( P-ip can be calculated for g 2 C using Wilson's operator product expansion (OPE) |TJj] 
for quark and gluonic degrees of freedom || (for in-medium calculations see e.g. |l|[n|). In 



the following we shall call the result of that calculation i? OPE . A second representation in 
the space-like region which has to match R OFE can be obtained from i? HAD by utilizing a 
subtracted dispersion relation. We find 

* + *_£Lls*^£> (2 . 6) 



Q 2 t { + Q 

with Q 2 := — q 2 S> and some subtraction constants Cj. 



Eq. ( |2.6|) connects hadronic with quark-gluon based expressions. In principle, for a given 
hadronic parametrization of i? HAD with free parameters this equation could be used to ex- 
tract information about these parameters. This, however, would require the knowledge of 
_R HAD (s) for arbitrary large s. In practice, the situation is such that one has a parametriza- 
tion for the current-current correlator for the energy region of the lowest hadronic resonance, 
but one usually has no model which remains valid for arbitrary high energies. In the dis- 
persion integral of ( |2.6| ) higher lying resonances are suppressed, but only by a factor 1/s 2 . 
Clearly, it is desirable to achieve a larger suppression of the part of the hadronic spectral 
distribution on which one has less access. To this aim, a Borel transformation P,|3"2]] can be 
applied to ( |2.6|) . For an arbitrary function f(Q 2 ) the Borel transformation is defined as 

f(Q 2 ) A f{M 2 ) (2.7) 

with 

Q2/iV=:M2=flxed V ' 

where M is the so-called Borel mass. Applying the Borel transformation to (|2.6| ) we finally 

get M 
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o+ 

We observe that higher resonance states are now exponentially suppressed. Note that the 
subtraction constant C2 of ( |2.6| ) has dropped out. The other one, ci, vanishes in vacuum. 
In a nuclear medium for a meson at rest, we incorporate the Landau damping term in the 
subtraction constant C\. This term comes from the absorption of a space-like meson by an 
on-shell nucleon. Having incorporated this term in C\ we avoid double counting by restricting 



the integration in ( |2.9| ) to the time-like region. For a detailed discussion of that point cf. 
One gets in the linear density approximation 

c 1 = /^. (2.10) 

Eq. ( |2.9|) is the QCD sum rule which we will utilize in the following. 

Having achieved a reasonable suppression of the energy region above the lowest lying 
resonance the integral in Q2.9|) is no longer sensitive to the details of the hadronic spectral 
distribution in that region. For high energies the quark structure of the current-current 
correlator is resolved. QCD perturbation theory becomes applicable yielding 

InLR HAD (s) = (\ + —) for large s. (2.11) 

071 \ 71 J 

These considerations suggest the ansatz 

IhlR had (s) = 0(s o - s) lmR RES (s) + &(s - s ) ( 1 + — ) (2.12) 

where sq denotes the threshold between the low energy region described by a spectral func- 
tion for the lowest lying resonance, Imi? RES , and the high energy region described by a 
continuum calculated from perturbative QCD. In the following we use a s (lGeV) ~ 0.36. 
Of course, the high energy behavior given in (|2.12| ) is only an approximation on the true 



spectral distribution for the current- current correlator. Also the rapid cross-over in ([2.1 2|) 



from the resonance to the continuum region is not realistic. However, exactly here the sup- 
pression factors discussed above should become effective making a more detailed description 
of the cross-over and the high energy region insignificant. The price we have to pay for 
the simple decomposition Q2.12D is the appearance of a new parameter so, the continuum 
threshold, which in general depends on the nuclear density. We will elaborate below on the 
determination of Sq. 



To study the content of (|2.9| ) for the vector and axial-vector channel we need the OPE 
for the l.h.s. In general, it is given by a Taylor expansion in 1/M 2 : 

£ope (m2)=E _!__ (213) 
In the following we present the formulae for the case of finite nuclear density p^. The 



vacuum case is easily obtained by p^ — > 0. For the vector channel one gets [0,0 (for 
details see also [E3 and references therein) 
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c V o = 


1 

1 

8tt 2 


/ rv \ 


c\ = 


0. 




<% = 


1/ 
24 \ 


7T / 




5 

~24 


m 3 N A 4 p N 



-m N 

56 
" 81' 



(2.14a) 
(2.14b) 
(2.14c) 

(2.14d) 



while for the axial-vector sector one obtains [@,[nj 

cf = c( for i = 0, 1, 2 , (2.15a) 
4 = ~^™%A 4 p N + ^a s (Of) . (2.15b) 



We neglect (unknown) condensates with dimension higher than 6 and some less important 
twist-4 condensates and ot s corrections (cf. [[0|[|l],|4[] ) . We also neglect perturbative contri- 
butions proportional to the square of the current quark masses. Note that all expectation 
values have to be taken with respect to the medium. We work here in the linear density 
approximation: 

(0)^(OU c + ^(N\0\N). (2.16) 

A single nucleon state is denoted by \N). It is normalized according to 

(N(k)\N(k')) = (2irf2E k 5(k-k'). (2.17) 

We defer the calculation of the in-medium expectation values of the scalar operators to 
Sec. |VT] and only discuss their vacuum expectation values here. For the gluon condensate we 
use a canonical value of M (—G 2 ) = (330 MeV) 4 . As compared to the gluon-condensate 

\ 7r / vac 

the influence of the two-quark condensate 

m q {qq) vac = --/X ( 2 - 18 ) 

is rather small (and is further diminished in a nuclear environment). Here f n = 93 MeV 
denotes the pion decay constant and m, the pion mass. While the values for gluon and two- 
quark condensate are fairly well known the knowledge about the four-quark condensates^ 

= ^ (feA^ - ^ 75 W) 2 ) + ^ ((«7,^ + E (2.19) 



and 



3 Note that the definition of (0^ A ) is chosen such that the factorization assumption would imply 
P\ IA ) « <W> 2 . 
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(Of) 



352 



(u^X a u - d^X a d) 2 



176 



(u^u + d^d) ^VW)- (2.20) 



lj)=U,d,! 



is very limited. Traditionally factorization is assumed which, however, probably underesti- 
mates its value. In the following we will use two values for the four-quark condensates to 
explore the sensitivity of the results: 



(C^4 )vac j {@4 )\ 



-292 Me V) 6 , (-281MeV) f 



[2.21] 



The larger value is chosen as to obtain an optimal agreement between QCD sum rule predic- 
tion and experiment for the p meson properties in vacuum. Finally the terms proportional 
to A 2 and A 4 in ( |2.14|,|2.15 ) stem from twist-2 condensates. They are obtained from the 
moments of the quark distributions in a nucleon 



0.9, A 4 = 0.12. 



We use A 2 

So far we have not specified for which values of M 2 we regard the sum rule Q2.9 ) to 
be valid. Note that in practice ( |2.13| ) is a truncated series in 1/M 2 . Clearly, if M 2 is too 
small the 1/M 2 expansion in ( [2.13 ) breaks down. On the other hand, however, if M 2 is 
too large the exponential in fl2.9| ) does not sufficiently suppress the intermediate and high 
energy part of IhlR had (s) given in ( |2.12| ). As mentioned above this suppression is important 
since the modeling of the region around the threshold so is rather crude. If these qualitative 
considerations are put on a more quantitative level one can define a so called Borel window 
for the masses M 2 in which the sum rule is valid (cf. e.g. [23]). Following |30| we determine 
the minimal Borel mass such that the last accessible contribution to the OPE ( |2.13| ), i.e. here 
the 1/M 6 term, amounts to 10% of the total OPE result: 



c 3 



ML 



0.1i? OPE (M^ m ) 



(2.22) 



The maximal Borel mass is chosen such that the continuum contribution to the r.h.s. of 
( p.9| ) does not become larger than the contribution from the resonance which we want to 
study, i.e. 



oo 

/^( 1+ ?) 9 < s 



so) e 



-8/M, 



2 

max 



>3C 

Jds IhlR res (s) 6(s - s) e 



-s/M, 



2 

max 



(2.23) 



As a guideline one can expect that M^ ax scales with the point where the average strength 
of Imi? RES (s) is located, i.e. with the resonance mass (squared). Hence for large (small) 
resonance masses the value of M^ ax will be large (small). It might appear that in some 
cases the Borel window between M 2 -,„ and M 2 ^ is rather small or even closed. Then the 
sum rule is meaningless. In practice the determination of the Borel window provides a 
quality check for the sum rule. 



III. FINITE ENERGY SUM RULES THE QUALITATIVE PICTURE 



The Borel sum rule ( |2.9|) is not the only sum rule which is used to connect hadronic 
and QCD based information. Inserting ( ^.12j ) in ( |2.9j ) and expanding the r.h.s. in powers 



10 



of 1/M 2 one can compare the coefficients of this expansion with the respective ones in the 
series on the l.h.s. given by ( [2. 13 ). This yields the finite energy sum rulesQ (presented here 



for the p meson for the vacuum case) 



1 S ° 

- fdslmR^ ES (s) - c^s = , (3.1a) 

71 J 


- Ids s Imi$ ES ( S ) + c y $ = c\ , (3.1b) 

7T J 2 


1 S ° 3 

- Ids s 2 Imi$ ES ( S ) - c v Q $ = 2c\ (3.1c) 

7T J 3 




where the coefficients of the OPE are given in fl2.14| ), evaluated in this section for p^ = 



The first two of these sum rules are utilized e.g. in jEJ. Obviously the expansion of 
the r.h.s. of ( |2.9| ) relies on the assumption that the Borel sum rule obtained by the simple 
decomposition ( |2.12|) is valid for arbitrary high values of M 2 . As pointed out above this is 



doubtful due to the limited knowledge of IhlR had (s) in the threshold region. Actually the 
sensitivity of the respective finite energy sum rule on the details of Imi? HAD (s) around sq is 
increasing when going from fl3.1a|) to (|3.1cj) . Thus it might be safe to extract information 
from the lowest finite energy sum rule(s). Utilizing higher ones, however, becomes more and 
more doubtful. This is the reason why we prefer to use the Borel sum rule. In addition, 
for the latter one a consistency check on its validity is provided by the determination of the 
Borel window.^ 

Nonetheless, the finite energy sum rules can be used to get a qualitative picture about the 
connection of the OPE side to resonance parameters like mass and width. In fact, Imi?y ES 
is a mass distribution. Therefore it appears natural to define the first two moments of this 
distribution, i.e. an average mass and a width via 

JdsslmR^ ES (s) 

m 2 := ^ (3.2) 

fds ImF$ ES (s) 
o 

and 



4 This derivation is actually oversimplified since it neglects the running of the coupling constant. 
For a rigorous derivation cf. [^] and references therein. 

5 A consistency check for finite energy sum rules might be obtained in the following way: Clearly 
the discontinuity between the resonance and the continuum region in ( |2.12| ) is unrealistic. It is only 
used to avoid new additional parameters. Introducing instead a smooth cross-over one can test 
the sensitivity of the finite energy sum rules on these new parameters which model the cross-over 
region. If a finite energy sum rule appears to be fairly insensitive to these new parameters it might 
be regarded as useful. 
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fds (s - m 2 ) 2 Imi?^ ES (s) 

a 2 m 2 := «— ^ . (3.3) 

JdslmR^ ES (s) 
o 

Obviously, the finite energy sum rules ( p.l|) can be used to connect these moments with the 
condensates (and the continuum threshold): 

m = — , (3.4) 

/ CqSq 



° l = — [^ + — -m* . (3.5) 
m 2 \ 3 c s J 

We can learn two things from these simple relations: First, the average mass is determined 
by the dimension-4 (gluon and two-quark) condensates and the continuum threshold while 
the dimension-6 condensates (here the four-quark condensate) influence only the width. 
Second, we do not have enough information at hand to determine all the phenomenological 
parameters. In our case at hand we have three of them, namely the continuum threshold 
So and the two moments rh 2 and a 2 . On the other hand we only have two equations for 
these parameters. Traditionally, the use of QCD sum rules is accompanied by an additional 
assumption, namely that the width is negligible. In this case the mass can be determined. In 
general, however, the best we can hope to gain are correlations between the free parameters. 
Being especially interested in mass and width we can vary sq and determine the correspond- 
ing values for rh and a. The result is shown in Fig. |I[ (For the four-quark condensate ( |2.21| ) 
we have chosen the larger value.) The most important thing to note here is that the width 
grows with rising mass. We will find correlations of this kind again and again throughout 
this work. The qualitative understanding of this correlation is obtained from the simple 
relations ( |3.4| , |3.5| ). We are reluctant, however, to draw any quantitative conclusions from 
the previous considerations. In principle we are interested in the properties of the vector 
and axial-vector resonances, e.g. their masses and widths as defined via Breit-Wigner type 
parametrizations. In general, these masses and widths are not identical to the moments rh 2 
and a 2 of the distribution IhlR res . In addition, as outlined above we doubt the quantitative 
reliability of the finite energy sum rules ( |3.1| ) due to their higher sensitivity to the details of 
the high energy behavior. Hence we prefer the use of the Borel sum rule (|2.9|) . For a very 



elaborate use of combinations of finite energy sum rules we refer to E2] 



IV. BREIT-WIGNER PARAMETRIZATION OF THE CURRENT-CURRENT 

CORRELATOR 

The only remaining question is how to parametrize Imi? RES (s) which enters the sum rule 
(2.9) via ( |2.12|) . Concerning the vector channel, experiments which determine e.g. e + — > 



7r + 7T suggest the parametrization 



lmR* ES (s)=irF p ^ (4.1) 
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Here F p determines the absolute height of the spectrum and S p denotes the spectral function 
of the p meson which we will specify further below. Concerning the axial-vector channel not 
only the a\ but also the pion shows up there. Hence the parametrization has to be extended 
to 

Imi?* ES (s) = «Ki ^ + nfj(s - ml) . (4.2) 

s 

The spectral functions are given by 

S ^ = L ( ^Ji S) tv m ■ ^ 

7T [s — m z ) z + s 1 (s) z 

Here m is the mass of the respective meson and T its width. We stress again that these 
Breit-Wigner parameters are not identical to the moments introduced in the last section; 
there is only a qualitative correspondence. We denote the on-shell width by 



7 = r(m 2 ) . (4.4) 



We have to use an s-dependent width in ( f4.3|) for the following reason: In the following we 
will vary 7 (and other parameters) over large ranges. As outlined above the sum rule ( [2.9| ) 
is insensitive to the modeling of the high energy behavior of Imi? RES (s). In turn, there is 
a high sensitivity to the low-energy part. Therefore, especially for large widths we have 
to make sure that at threshold the spectral function shows the correct behavior. On the 
other hand, we don't want to overweight our parametrizations with too many independent 
parameters. Hence, we are aiming at simple parametrizations which reproduce correctly the 
threshold behavior. 

In the vacuum the width of the p meson is governed by the decay into two pions. We 
use the following parametrization: 



rr ay (s) = 7p^KfFk ( 4 - 5 ) 




with the momentum of the pions in the rest frame of the decaying p with invariant mass 

p:: i (s) = (s-4ml)^/2. (4.6) 

Concerning the a\ meson in vacuum its width is dominated by the decay into rho plus pion. 
For simplicity we neglect the width of the rho meson here and use 



rtr y (s) = lai^^4-Ms - (m p + m n ) 2 ) (4.7) 



with the momentum of pion and rho in the rest frame of the decaying a\ with invariant mass 

pSOO = Us - K + ^) 2 )(* - K - ^) 2 )1 1/2 /(2V5) • (4.8) 
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In a nuclear environment a presumably rather sizable collisional width (cf. |39] and references 
therein) has to be added to the decay width. The lowest threshold for a (axial-)vector meson 
collision with a nucleon is given by the formation of pion plus nucleon. We assume that the 
threshold behavior is dominated by the lowest accessible partial wave. For the p meson this 
is an s-wave leading to 

rf\s) = % f /~f/0 B{* - O • (4.9) 



For the a\ meson it is a p-wave: 



rffOO = Toi 4- 1 r 7/°2 ) B(s-ml). (4.10) 
mix V-mi/m z al J 



1 — m 2 I s 



3/2 



For the vacuum case r p / ol (s) in Q4.3| ) is given by ( Op and (|4.7|), respectively. For the case 
of nuclear medium we restrict ourselves to the two extreme possibilities that the width is 
either dominated by decays or by collisions. Hence we explore the two cases that r p / al (s) is 
either given by (Pfl/(pfl or by flPQ/flETq). 

We will treat F, m, 7 and also the continuum threshold Sq (cf. (|2.12|) ) as free parameters. 
The aim is to find out how the sum rule ( |2.9| ) constrains these parameters. As a general 
rule one can at best determine as many parameters of the hadronic side of the sum rule as 



one has powers in 1/M 2 on the OPE side For the latter we have given in ( |2.14|j2TT5|) 
four orders in powers of 1/M 2 . However, the perturbative part ((1/M 2 ) part) has already 
been used to determine the high energy behavior in Therefore at best only three 



parameters of the hadronic spectral distribution can be determined from the sum rule ( |2.9|) .f| 
On the other hand, we have for each meson four free parameters in the parametrization 
(CT,P1,P1,PD. In the traditional sum rule approach 0jg|EJ|3gg| the width of the 



respective meson resonance is neglected (narrow width approximation). In this case the 
number of free parameters reduces to three and the sum rule gains predictive power. This 
however means that in addition to the QCD input represented by the OPE one needs further 
knowledge to extract predictions from QCD sum rules. In vacuum, this additional input 
comes from experiments which tell us that e.g. the p meson indeed is a well-defined resonance 
with a width considerably smaller than the mass. In contrast, for the in- medium case it 
is so far not clear if the pronounced peak structure of the p meson survives in a nuclear 
surrounding or if it is washed out [ [27|j35H 39l e.g. by its coupling to resonance-hole states. 
The ax meson already has a large vacuum decay width. In addition, its mass is so high that 
it is hard to achieve a clear separation between Imi? RES and the high energy continuum in 
( |2.12j ) 0. We will come back to that point below when discussing the results of our QCD 



sum rule analysis for the a\. Like for the p meson, medium effects in addition presumably 
lead to an additional broadening of the a\. To study the influence of the widths of the vector 



6 This is also true if the finite energy sum rules were used instead of the Borel sum rule. The 
three sum rules (3.1) provide three constraints on the hadronic spectral distribution. See also the 
discussion in Sec. ITTl. 
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and axial-vector meson on the results extracted from QCD sum rules we will refrain from 
neglecting 7 and proceed with our general parametrization (|4.3| ). 



V. RESULTS FOR VACUUM 

We shall now explore which values of mass and width of the p/a± meson are compatible 
with the sum rule ( |2.9| ). For that purpose we vary the values for m and 7 in large ranges. 
For a given pair of mass and width we tune the remaining parameters F and Sq such that the 
agreement between left and right hand side of (|2.9| ) is best. The resulting minimal deviation 
d between l.h.s. and r.h.s. is a measure for the compatibility of the chosen mass-width pair 
with the sum rule, i.e. if d is sufficiently small one might conclude that the chosen pair of mass 
and width is allowed by QCD sum rules. We regard the sum rule to be approximately valid 
in the range given by the Borel window introduced above. Hence we define the deviation d 
as an average over this window (see [f^J for further details). 

For the p meson Figs. § and |3| show the allowed ranges for mass and width for the vacuum 
case for the two different values of the four-quark condensate given in ( |2.21| ). Obviously, there 
is not only one point where the sum rules is reasonably fulfilled but a whole band of allowed 
mass-width pairs. Figs. [2] and ^| qualitatively resemble Fig. [TJ: the band of allowed mass- 
width pairs describes a correlation where the width rises with rising mass. Quantitatively, 
however, the differences between Fig. [I] and Figs. [|, ||| are large stressing again that the 
distribution moments defined in Sec. [TTT] are not identical to mass and on-shell width of the 
Breit-Wigner type spectral functions ( |4.3|) . We have also included the experimental point for 
the p meson in Figs. |2| and [| Obviously the value of (—292 MeV) 6 provides an optimal choice 
for the four-quark condensate. A smaller value for that condensate shifts the band of allowed 
mass- width pairs either upwards or to the left (or both). From our qualitative considerations 
of Sec. m we expect that the four-quark condensate mainly influences the width (cf. ( |3.5|) ). 



Hence a smaller value of the four-quark condensate is supposed to increase the width,[]i.e. to 
shift the band upwards. We will further clarify that point when discussing the results for the 
a\ meson. Note that also for the smaller value of the four-quark condensate of (—281 MeV) 6 
the deviation d is still reasonably small for the experimental point (0.46%). As discussed 
above the results of Figs. |2| and |] are only meaningful if the Borel window is reasonably 
large. For our choice for the condensate values leading to Fig. |2] we get M^ in = 0.71 GeV 2 . 
The value for M^ ax depends on the resonance parameters. As already mentioned small 
values for the resonance mass lead to small values for M^ ax . For the case at hand we find 
M^ ax > 1.5 GeV 2 for all mass-width pairs lying in the inner band shown in Fig. [| and to the 
right of it. We regard that as a reasonably large Borel window. For Fig. |3| the corresponding 
values are M 2 in = 0.65 GeV 2 and M 2 ax > 1.4 GeV 2 . 

Concerning the a\ meson the corresponding mass- width correlations are shown in Figs. |] 
and [|. Qualitatively we find again the same correlation between masses and widths. How- 
ever, a tendency is visible that the sum rule supports large values of the width. We will find 



7 Note that the four-quark condensate enters the sum rule for the p meson with a negative sign, 
cf. (pl4d| ). 
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that this tendency increases for the in-medium case. Comparing Figs. ^ and |5| we find that a 
decrease in the four-quark condensate shifts the band to some extent to the left but merely 
downwards. To understand that finding we note that the four-quark condensate enters with 
a different sign in the two sum rules for p and a%, respectively; cf. fl2.14d| , [2~.15b|) . Therefore 
we expect an upward shift with decreasing four-quark condensate for the p meson, as dis- 
cussed above, and a downward shift for the a%. Figs. f| and || support our considerations. 
Comparing the results with the experimental values for mass and width of the a\ meson we 
find that the smaller value for the four-quark condensate (Fig. [5]) appears to be much better 
suited — quite opposite to the case of the p meson where the larger value provides a better 
fit. In principle, there is no fundamental reason why the four-quark condensates for p and 
ai should be exactly the same. Only if the factorization assumption strictly holds the two 
quantities defined in ( |2.19| ) and ( |2.20| ) coincide. Still, however, the values (—292 MeV) 6 and 
(—281 MeV) 6 are rather close. For the succeeding in-medium calculations we will choose the 
respective better value, i.e. (—292 MeV) 6 for the p meson and (—281 MeV) 6 for the a%. We 
note that for the a\ vacuum case a further reduction of the four-quark condensate does not 
improve the agreement between the sum rule and the experimental results as the band is 
further shifted downwards and not to the left. Hence the agreement between sum rule and 
experiment appears to be better for the p than for the a\ meson. This rather old finding |2j 
is most likely due to the fact that the separation between the resonance and the continuum 
region is better realized in the vector channel. If the resonance appears to be closer to 
the continuum the sum rule is more sensitive to the details of the modeling of the nearby 
transition to the continuum. Such details are necessarily rather crude in our ansatz Q2.12|) . 
Finally we present the results for the Borel window for the preferable parameter choice 
(Fig. |): We find M^ in = 0.71 GeV 2 while all mass-width pairs enclosed by the dashed line 
obey M^ ax > 1.4 GeV 2 . We note that we have not included the pion branch in the determi- 
nation ( |2.23|) of the maximal Borel mass to make sure that we really learn something about 
the properties of the a\ meson. 



VI. IN-MEDIUM CORRELATORS AND p-A 1 MIXING 

Next we turn to the case of finite nuclear density. As pointed out by several groups 
[e.g. [ f4^[45|| ) the interaction of the p meson with the pion cloud of the nucleons induces a 



mixing of the p with its chiral partner the a\ meson. This means that e.g. the (possibly 
medium- modified) ai-peak shows up in the spectral distribution lmRy AB of the vector corre- 
lator and vice versa. Suppose now that one would ignore that multi-peak structure and still 
parametrize lmR v AD with only one peak according to ( |2.12jJ4.1j ). In a QCD sum rule analysis 
one has only access on certain mass averages of the spectral distribution on account of Q2.9|). 
Hence with a one-peak structure ansatz one would translate certain in-medium changes of 
the OPE side to changes of mass and width of this peak which in reality, however, are caused 
by the appearance of other distinct peaks. This would be rather misleading. Indeed, for 
the comparable case of finite temperature T it has been shown [|[ that the masses of p and 
Oi, if understood as the positions of peaks in spectral distributions, do not change in the 
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linear (pion) density approximation.^ Only if the notion of mass is used with a different 
meaning (e.g. in the spirit of Sec. |T| as the first moment of a spectral distribution) it would 
be correct to attribute an in-medium 0(T 2 ) mass shift to this "mass". For considerations 
beyond the 0(T 2 ) approximation we refer to |16| , |T7 . 



Concerning the present work the mass m of a resonance (which shows up in the spectral 
distribution of the current correlator) is defined via (|4.3|) . For small width it gives the peak 
position of the resonance. In the following for our case of finite nucleon density we also try to 
account for the multi-peak structures caused by mixing of vector and axial- vector currents. 
If we only used the sum rule ( ^.9[ ) and introduced more than one peak e.g. in ( |4.1| ) we would 
have too many free parameters to draw any meaningful conclusion. Therefore (as for the 
case of finite temperature J|) the key idea is to isolate the contribution of the respective new 
in-medium peak(s) also for the OPE side. In this way one obtains sum rules for non-mixed 
correlators (in the following called "bare" ) which can be analyzed with the one-peak ansatz 
(plus continuum, of course). These (in general unobservable) bare correlators mix to yield 
finally the "full" in-medium correlators. The imaginary part of the latter in principle can be 
observed (cf. e.g. Q2.5|)). We note that for the case of finite temperature the bare correlators 
coincide with the vacuum correlators in the linear density approximation |4]]. As we shall 
see in the following, things are not so simple for the case of finite nucleon density. 

To account for the interaction of the nuclear pions with the vector- and axial-vector 
currents nuclear matter is separated into a Fermi gas of bare nucleons plus soft pions, 
schematically |44 



|tf > = * A \A) + E <\ A ^) + E <M*a* b ) + ... (6.1) 

a. a,b 

where |^) denotes the full nuclear matter state vector while \A) denotes the bare one. The 
current- current correlator ( p.l|) evaluated with respect to |\?) can now be decomposed into 



a bare correlator, i.e. a correlator with respect to \A), and a part involving the interaction 
with (soft) pions: 

IV = nj, + n; u . (6.2) 



The latter one is approximately evaluated using soft pion techniques (cf. also [151) and taking 
into account up to two pions in the initial and/or final state. One gets 



(A7T a \Tj,(x)j u (0)\A7r b ) ~ (An a 7r b \T Jtl (x)U0)\A) ^ (A\Tj,(x)U0)\An a n b ) 

-1 

j - 



\a\ \Ql\Q b 5 ,Tj,(x) 3u (0)}} \A), (6.3) 



with the isovector axial charge 



(X; = /r/ : Vr(.r)~„-^r(.r). (6.4) 



8 Note that an 0(T 2 ) modification at finite temperature corresponds exactly to 0(pn /tun) at 
finite baryon density (cf. e.g. Eq. (3) in Q). 
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To calculate the commutators in (|6.3| ) with the currents (|2.2| ) and ( |2.3|) it is useful to gen- 
eralize the latter to the full isospin multiplets 



(6.5) 



In fact, on account of 32 



ie abcyc 



(6.6) 



the vector and axial-vector currents are mixed by their interaction with the nuclear pions. 
Finally the expectation values in (|6.3| ) have to be weighted by the density of pions in the 
nuclear medium. For the dimensionless quantities defined in ( |2.4| ) one ends up with (see ||4 
for details) 



Rv(g 2 ) = rW) - aRW) - Ra{Q 2 )) 

R A ( q 2 ) = R b A (q 2 ) - aRW) - RW)) 



(6.7a) 
(6.7b) 



where Ry/ A denotes the respective correlator with respect to a system of bare nucleons, 
i.e. without their pionic cloud. The mixing parameter is given by 



3 f^ml 



(6.8) 



Here we have introduced 



(6.9) 



where N w denotes the scalar number of pions in the cloud surrounding the nucleon. a 1 ^ 
contributes to the nucleon sigma term given by [p4j 



m N 



(N\qq\N). 



(6.10) 



Having split up the nucleons into bare nucleons plus a cloud of soft pions we have to disen- 
tangle the nucleon sigma term correspondingly ||45|| : 



'N 



(6.11) 



At present, the value and even the sign for is not a settled issue. We take a positive 
value of 25 MeV. This is in agreement with but in contrast to [11] where a negative 
value has been used. Our choice for ajj is motivated by the fact that this ensures that the 
two correlators Ry and Ra become degenerate at some finite density: 



Ri 



R A = (l-20(R b v -R b A ). 



(6.12) 



A negative value for ajf would lead to anti-mixing, i.e. in this case the nuclear pions would 
work against chiral symmetry restoration. 
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Next we will perform a sum rule analysis for the correlators with respect to the system 
of bare nucleons Ryi& at nuclear saturation density = 0.17/fm 3 . Concerning the vector 
channel this is the essential difference as compared to our previous work |23[ where we have 
analyzed the sum rule for the full in-medium correlator Ry. Note that Q2.14Q , (|2.15| ) are 



valid for both full and bare correlators. The difference appears in ( |2.16| ) in the calculation 
of the expectation value with respect to bare or full nucleons, respectively. In practice, the 
difference manifests itself only in the different handling of the four-quark condensate. All 
OPE contributions in (|2.14| ) and ( |2.15| ) except from the two- and four-quark condensates 



come from chiral singlet operators. They do not distinguish between bare nucleons and 
nucleons dressed by soft pions. Hence their evaluation is standard [P^fl . Things are different 
for the two- and four-quark condensates. Concerning the two-quark condensate, its in- 
medium change (in linear density approximation) is determined by the nucleon sigma term 
( |6.10| ). To calculate the in- medium change with respect to bare nucleons we have to take 
into account only a b N as defined via ( |6.11| ): 



(qq)b 1 QnPn 

-i-^T-T- ( 6 - 13 ) 



(qq) vac f%m% 

In practice, the value of m g (qq) is rather small (as compared to the gluon condensate) and 
further diminished in the medium. In contrast, the four-quark condensate is numerically 
important. In lack of a better access to the value of 0\ at finite density we made in [^] the 
common assumption that this four-quark condensate scales with the density like the square 
of the two-quark condensate 

In the comparable case of finite temperature (i.e. in a hot pion gas), however, it was shown 
46]] that such an assumption is wrong. This can be traced back to the fact that in the 



in 



presence of pions the two-quark condensate behaves different as compared to the four-quark 
condensate due to its different transformation properties with respect to chiral transfor- 
mations. This suggests that also at finite nucleon density the scaling assumption fl6.14|) is 



doubtful due to the presence of virtual pions. In the present work we have explicitly taken 
into account the contribution from the pion cloud of the nucleons. In this way we have 
expressed the full correlator in terms of the bare correlators. We now assume the scaling 
property (|6.14j) only for the condensates with respect to a system of bare nucleons. It takes 
the form 



(6.15) 



where for consistency we have to take the bare nucleon sigma term a b N = — ajj ~ 20 MeV 
instead of the full one <jn ~ 45 MeV. In the following we use this scaling assumption (|6.15|) 



for both (OX)b and (Of) b . 

In view of the uncertainty connected with the four-quark condensate it clearly would 
be fortunate to use sum rules which do not involve it. Indeed, in the first two finite 
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energy sum rules ( |3.1a|J3.1fj|) were used. However, with the same parameter set (m p , 7, F, sq) 



characterizing the hadronic correlator the two sum rules are capable to determine only two 
of these four parameters. To further restrict the parameter space additional information 
is required. In [|3|] it is suggested that the threshold s is connected to the scale of chiral 



symmetry breaking. Thus, the choice is either to make assumptions about the four-quark 
condensate or about the threshold parameter. As outlined above we prefer to work with the 
Borel sum rule instead of the finite energy sum rules due to the larger sensitivity of the latter 
to the high energy modeling. In this case we cannot get rid of the four-quark condensate. 

The sum rule analysis proceeds along the same lines as described above for the vacuum 
case. We analyze the Borel sum rules for the p as well as the a\ meson placed in a medium 
of bare nucleons. For the time-like part of the correlator in the vector channel we use again 
a single resonance parametrization of type ( |4.1| ). For the axial- vector channel we recall that 
there is a pion branch in addition to the a\. Here there is an additional change in the 
medium due to a change of the pion decay constant in nuclear matter. We replace f% in 
O by 

L ~ /ttT-t — (6-16) 

(Wvac 

where we have utilized the in-medium version of the Gell-Mann-Oakes-Renner relation 
( p,18| ) and neglected a possible in-medium change of the pion mass. By using ( }4.2| ) with 
the replacement (|6.16|) we also assume that the delta-type spectral function of the pion is 



not significantly smeared out. In fact, pion properties are expected to change drastically in 
nuclear matter due to the strong coupling to Delta-hole states [f47|. However, this p-wave 
coupling is not important here since we deal with correlators which are at rest with respect 
to the nuclear environment. Working in the linear density approximation, i.e. neglecting 
Fermi motion, the pions do not couple to the nucleons by exciting Deltas. In the context 
of the (necessary) approximations involved it seems reasonable to work with a pion spectral 
function which is neither shifted nor broadened. 



VII. RESULTS FOR NUCLEAR MATTER 

The results of the sum rule analysis are shown in Figs. || [7|, and |8|. As for the vacuum 
case we find bands of allowed mass-width pairs. The interesting point is how these bands 
have changed as compared to the respective vacuum case. As already mentioned concerning 
the choice for the vacuum four-quark condensate we have used the respective better value, 
i.e. (—292 Me V) 6 for the p meson and (—281 Me V) 6 for the a\. Hence we have to compare 
Figs. |, to Fig. I and Fig. || to Fig. [|. 

For the p meson we have explored two possibilities to parametrize the energy dependence 
of the width. The same parametrization (|4.5| ) as for the vacuum case, i.e. with the two-pion 
threshold, is used to obtain the results depicted in Fig. Here the mass-width band is 
shifted to the left as compared to the vacuum case (Fig. 0). Obviously, the in- medium 
change of the condensates calls for more strength at lower invariant masses as compared 
to the situation in vacuum. This can be accomplished either by a lower peak mass or by 
a larger width. For Fig. [7] the parametrization ( |4.9| ) with the one-pion threshold is used. 
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For very small width the masses allowed by the sum rule analysis agree for Figs. |6| and 
|7]. This is of course due to the fact that the details of the off-shell parametrization of the 
width do not matter if the width is sufficiently small. For larger widths, however, there 
appear large differences between Figs. |] and |7|. The band in Fig. [7| is much less steep. This 
is easy to understand: The demand for more strength at lower invariant masses is easier 
to fulfill if there is already some strength below the two-pion threshold. Therefore, for the 
same peak mass the on-shell width can be smaller if the one-pion threshold is used instead 
of the two-pion threshold. The large differences between Figs. ^| and |7| stress again that the 
Borel sum rule is very sensitive to the low-energy behavior of the spectral distribution. More 
generally one has to realize that peak mass and on-shell width are sufficient to characterize 
a spectral function only if the width is not too large. For large width details of the spectral 
shape become important, in the case at hand especially the details in the low-energy region. 
Concerning the Borel mass window we get M^ in = 0.64 GeV 2 . For all mass- width pairs in 
the inner bands of Figs. ^| and |7] and to the right of it we find M^ ax > 1.2 GeV 2 . It is a 
generic finding of sum rule analyses that the Borel window shrinks to some extent when 
changing from the vacuum to the in- medium case |[23|| . Still we regard the Borel window to 
be large enough to draw conclusions from the analysis. 

Turning to the ai meson Fig. |^ shows that the tendency of the sum rule to support large 
values of the width (cf. Fig. [5]) increases in medium. (Note that the width-scale in Fig. |8| 
differs from the previous figures.) The minimal Borel mass is given by M^ in = 0.61 GeV 2 
while M^ ax > 2.0 GeV 2 for the whole relevant mass- width range. For the width we have used 
parametrization ( |4.10| ) with the one-pion threshold to obtain Fig. || We have also analyzed 
the sum rule using instead ( £7 ) with the rho-pion threshold. In this case we did not find 
any mass-width pair with a deviation d less than 1%. Therefore, we do not show a plot for 
the latter case. In fact, even if we allowed for larger values of d there would be no sign for 
the desired mass-width band. We conclude that parametrization ( |4.7|) is incompatible with 
the in-medium sum rule for the a>\ meson. Obviously this sum rule demands for a spectral 
distribution which is smeared out over a large invariant mass range. Restricting the mass 
range by a (rather high) lower limit of m n + m p ~ 0.9 GeV appears to be insufficient to fulfill 
the sum rule. Instead, the one-pion threshold (caused by nucleon-ax collisions) does the job. 
Of course, also other scenarios which might fulfill the sum rule are conceivable. As we have 
learned from the previous in-medium sum rule analysis the spectral function of the p meson 
gains strength at lower invariant mass (either by a lower peak mass or by a larger width). 
Thus the effective threshold for the in-medium decay of the a\ meson into rho plus pion is 
probably lowered as compared to the vacuum case. Such a scenario might also be in line 
with the sum rule. We have not explored this possibility in further detail since we would 
have to introduce a couple of new free parameters to model e.g. the successive in-medium 
decays a± — > p + n — ► 3n. As already noticed when discussing the results for the p meson 
we find also for the a x that the QCD sum rule is sensitive to the threshold modeling of the 
width provided that the on-shell width is not too small. This fact has not been sufficiently 
taken into account in leading to results for the a\ meson which differ from the ones 
presented here. We also deduce from Fig. |8] that the nominal peak mass is shifted to higher 
values. It is important to note that this shift is not in contradiction to chiral symmetry. 
What we have studied so far is the behavior of the correlators under the influence of the 
system of bare nucleons. Chiral symmetry only demands that the full correlators of the 
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vector and axial-vector channel become degenerate at high enough density. This is achieved 
by the pions as expressed by the mixing formula ( |6.12j ) no matter how the a% mass changes 
under the influence of the bare nucleons. 

To visualize the effect of mixing we now turn to the full in-medium correlators obtained 
from the bare ones via ( |6.7| ). As already pointed out the sum rule analysis in general is not 
capable to pin down both in-medium mass and width for the respective meson. Nonetheless, 
to illustrate the effect of mixing we arbitrarily take one pair of values for the p meson from 
the inner band depicted in Fig. [7| and the optimal pair for the a\ meson. We choose: 

m p = 0.77 GeV, 7 p = 0.21GeV, s p Q = 1.15GeV 2 , F p = 1.1 • 1(T 2 GeV 4 ; (7.1a) 
m al = 1.86 GeV, j al = 1.67 GeV, s^ 1 = 3.84GeV 2 , F al = 0.36 GeV 4 . (7.1b) 

The results for the imaginary part of the full in-medium vector and axial-vector correlators 



are shown in Figs. ^ and |10] together with the respective contributions from the bare correla- 
tors. In the figures we have not included the delta function type contributions from the pion 
and from Landau damping. Note that the shoulder of the p contribution at low invariant 



mass is caused by the 1/s factor present in (3~1). Especially the axial- vector correlator shows 



that in the medium the resonance peaks are no longer higher than the high energy contin- 
uum. Recalling that the modeling of the onset of the continuum is rather crude one might 
further soften the crossover regions leaving basically no room for distinct peak structures. 
In total, we see a clear sign that the in-medium spectral distributions get washed out with 
increasing density. 



VIII. SUMMARY 

We have presented a QCD sum rule analysis for the in-medium current-current correlators 
with the quantum numbers of p and a\ mesons. For the medium we have chosen the case of 
a Fermi gas of nucleons at vanishing temperature. For comparison we have also presented 
the vacuum sum rule analysis for p and a\. As possible in- medium changes for the spectral 
distributions of both correlators we have allowed for mass-shifts, peak broadening and also 
mixing. The latter effect has not been considered in previous analyses |18| , ^9| , |23 1. 



In the QCD sum rule only a mass-averaged quantity involving the respective spectral 
distribution enters. In general the sum rule is not capable to pin down the full information 
present in the spectral distribution like the number of peaks, their positions, widths and 
heights. Therefore one needs an ansatz for the spectral distribution with some free parame- 
ters. These parameters can be constrained by the requirement that the spectral distribution 
has to fulfill the sum rule. The simplest ansatz consists of only one peak (in the low energy 
regime) and neglecting its width. In this case the sum rule for the p meson demands an 



in-medium mass shift towards smaller masses p^JT9|1 . At normal nuclear matter density the 



mass has to drop by roughly 16%. This, however, is a model dependent statement since 
a specific choice for the form of the spectral distribution has been made, namely only one 
peak with vanishing width. In fact, if also the width is included but still using a one-peak 
ansatz it is already impossible to fix both the mass and the width. One obtains a band of 
allowed mass- width pairs p3[ . At small width this band of course has to start at the mass 



determined in flT8"| , |I§[ — provided the same condensate values are used (cf. [T23[] for details). 
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For larger widths the allowed masses also increase. With increasing nucleon density the 
mass-width band is shifted to the left, i.e. for a given value of the width to smaller masses. 

If the true in-medium spectral distribution possesses more than one peak but in the 
ansatz for that quantity only one peak was present then certain in-medium changes caused 
by the additional peaks would be erroneously attributed to changes of mass and/or width of 
the single peak. Therefore, in the present work we have made a further step to consider all 
possible in-medium modifications. The mixing of the vector with the axial-vector correlator 
by the virtual nuclear pions has been included. Using soft pion techniques this mixing can 
be approximately calculated for the correlators irrespective of the choice for the correlator 
momentum q. Therefore the correlators can be decomposed into a superposition of "bare" 
correlators for both the OPE side (q 2 <C 0) and the spectral distribution (q 2 > 0). In this 
way we have split off the mixing phenomena (at least the one induced by nuclear pions) from 
the sum rule analysis. Successively the latter has been applied to the bare correlators. As 
compared to previous works we have obtained a less drastic in-medium shift of the allowed 
mass-width band for the p branch. For the a\ branch we have found that the sum rule 
is better fulfilled for not too small values of the width. The preference of large widths 
increases with increasing density. We have also seen that for both the p and a\ meson at 
larger widths the respective sum rule is sensitive to the threshold modeling of the energy 
dependence of the width. Therefore, it appears to be insufficient to characterize a peak in 
the spectral distribution by its position and on-shell width. A reasonable modeling of the 
energy dependence of the width is important to obtain meaningful results from a QCD sum 
rule analysis. 

As already mentioned we have found that there is a smaller in-medium shift of the mass- 
width band for the p meson as compared to the old analysis which did not include the 
mixing phenomena. At first glance this is a striking result. Suppose that the old analysis 
is polluted by an additional a\ peak, i.e. a peak at large invariant masses. First of all, one 
might argue that this additional peak is suppressed by the exponential function appearing 
in ( |2.9| ). Therefore this peak should not drastically modify the analysis for the p meson as 



generally assumed for all details of the high-energy part. Note that Figs. |9| and ^ show 
that in general the position of the a>i peak is beyond the onset of the continuum of the p 
branch, cf. also (|7. 1| ) . Let us ignore that point for a moment and face a second apparent 
contradiction of our results to naive expectations. Consider the influence of the additional a\ 
peak on the determined mass (for a given width) . Clearly the mass which we have attributed 
to the p meson in the old analysis would be a weighted average of the true p peak and the 
a\ peak. Therefore the true p peak should be at smaller mass values as compared to the 
result of the old analysis. In the present analysis, however, we have found just the opposite 
result. The solution to these apparent contradictions is the appearance of the pion in the 
axial- vector branch. It is not the high lying a\ but the low lying pion which has dominantly 
pollutes the old analysis. Since the pion is much lighter than the p the former has caused 
a too large shift of the mass-width band to the left in the old analysis. The a>i peak plays 
only a subdominant role as expected from modifications in the high energy regime. 

After the decomposition of the full correlators into bare ones we have tacitly assumed 
that for the latter the one-peak ansatz for (the low energy part of) the spectral distribution 
is reasonable. In fact, also this assumption appears to be questionable due to the coupling 
of the p-nucleon system to resonances (cf. e.g. PH|35| and references therein). Especially 
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the apparently sizable coupling to the D13 resonance iV*(1520) might create an additional 
peak at an invariant mass of roughly 580 MeV, i.e. below the vacuum p peakP] — provided 
that the resonance mass does not change in the nuclear medium. Therefore e.g. our result 
for the in-medium mass-width band of the p branch might still be polluted by the influence 
of additional distinct peaks. We have refrained from including more than one peak in the 
ansatz for the spectral distribution since we would have been forced to introduce much more 
free parameters which cannot be determined from the sum rule. Further work is necessary 
to separate the different peaks. On the other hand, we have already seen that the present 
sum rule analysis supports broad in-medium spectral functions. It may appear that the 
baryonic resonances also melt in a nuclear environment (cf. e.g. and references therein) 
and leave no sign of distinct peaks in the mesonic channels. 

Concerning the included mixing effect of the vector and axial- vector channel the amount 
of mixing is determined by the value for ajf. As we have already mentioned even the sign 
of that quantity is still subject to discussion. The influence of a 1 ^ is two-fold. First, it 
yields the amount of mixing on account of ( |6.7|) . Here a change of crjy does not influence the 
sum rule analysis for the bare correlators but it does influence the final result for the full 
correlators exemplified in Figs. |9] and [1^. Especially if cr]^ was negative one could no longer 
guarantee that the sign of the imaginary parts of the full correlators remains positive for all 
values of the invariant mass y/s. In this case the approximations which have led to ( |6.7|) 
should be revised. 

The second place where ajj appears is the in-medium dependence of the quark conden- 
sates. Most importantly, it enters the scaling assumption (|6.15| ) for the four-quark conden- 
sate via ( |6.11| ). We have made this assumption in lack of any better, more fundamental 
approach. One might corroborate our choice by the expectation that factorization roughly 
works as long as there is no fundamental (symmetry) principle which is in contradiction 
to factorization. In our case the interaction with the pions is determined by chiral sym- 
metry. Here factorization breaks down as proven for the case of finite temperature in [0rj 



Having split off the pions one might expect that the factorization assumption works for the 
interaction of the vector and axial-vector mesons with the bare nucleons. Nonetheless it 
is important to figure out how important the influence of the four-quark condensate and 
especially of its in-medium change actually is: At a typical value for the Borel mass of 
M « 1 GeV we find the following in-medium changes on the OPE side of the p meson sum 
rule at nuclear saturation density: The important changes are induced by the twist-two 
condensate with dimension four, 2.7 • 10~ 4 , and the in-medium change of the four-quark 
condensate, 1.5 • 10 -4 . All other changes are an order of magnitude smaller: —3.6 • 10~ 5 from 
the in-medium change of the gluon condensate, 1.3 • 10 -5 from the in- medium change of the 
two-quark condensate, —2.7- 10~ 5 from the twist-two condensate with dimension six. There- 
fore even if the change of the four-quark condensate was completely neglected roughly two 
thirds of the in- medium change would still persist. Thus, the results obtained here would not 



Note that only resonances with an s-wave coupling to N-p can contribute for the case at hand 
since we have considered p mesons at rest with respect to the nuclear medium. On account of the 
linear density approximation we have neglected the Fermi motion of the nucleons. 
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drastically change. The vacuum contribution of the four-quark condensate is —4.8 • 10" 4 . If 
the four-quark condensate vanished completely at nuclear saturation density the in-medium 
change would be twice as large as the one used in our analysis. One should realize that the 
previous considerations are somewhat oversimplified as the dependence on the Borel mass 
is different for the OPE contributions from the dimension-four condensates as compared to 
the ones with dimension six. In total, however, we do not expect qualitative changes of 
the picture presented here even if the in-medium behavior of the four-quark condensate was 
completely different. 

In spite of the large uncertainties connected with the QCD sum rule method at finite 
baryon density we regard it as a useful complement to purely hadronic approaches since it 
connects in a unique way hadronic with quark-gluon degrees of freedom. We hope that the 
mentioned uncertainties can be removed step by step in the future. 
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FIGURES 
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FIG. 1. Correlation between "mass" m and "width" a for the p meson as obtained from the 



finite energy sum rules (3.4,3.5). See main text for details. 
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FIG. 2. Deviation d as a function of width and mass of the p meson for vacuum. For the 
four-quark condensate a value of (— 292 MeV) 6 has been used. The full lines border the region 
of QCD sum rule allowed parameter pairs with d < 0.2%, the dashed lines border the allowed 
region for d < 0.5%. The cross marks the experimental values for mass and width of the p meson 
(including the (very small) error bars according to J50[|). 
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FIG. 4. Same as Fig. | for a\ meson in vacuum; value for four-quark condensate: (—292 MeV) . 
The cross marks the experimental values for mass and width of the a\ meson (including the error 
bars according to |5C 
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FIG. 5. Same as Fig. |2| for a\ meson in vacuum; value for four-quark condensate: (—281 MeV) 6 . 
The cross marks the experimental values for mass and width of the a\ meson (including the error 
bars according to ) . 
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FIG. 6. Same as Fig. || for p in system of bare nucleons at normal nuclear matter density. For 
the calculation of the width the two-pion threshold is adopted. See main text for details. 
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FIG. 8. Same as Fig. ^ for ai in system of bare nucleons at normal nuclear matter density. 
For the calculation of the width the one-pion threshold is adopted. See main text for details. Note 
the different scale for the width as compared to previous figures. 
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V~s [GeV] 

FIG. 9. Imaginary part of the full in-medium vector correlator (full line) as a function of 
y/s = \f(p- The contribution of the bare vector (axial- vector) correlator is the given by the dashed 
(dotted) line. See main text for details. 



Q 

< 



0.08 - 



0.06 - 



0.04 - 



0.02 - 







~ i i i i i i r 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

[GeV] 

FIG. 10. Imaginary part of the full in-medium axial-vector correlator (full line) as a function 
of yfs = \J~cp. The contribution of the bare axial-vector (vector) correlator is the given by the 
dashed (dotted) line. See main text for details. 
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